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Forward-backward correlation in pp collisions is studied in an approach that empha- 
sizes the partonic scattering angles and circumvents the intractable problem related to 
■ the transverse momenta that are low. Assuming the back-to-back scattering of partons 

to be the origin of hadronic correlation, the properties of forward-backward multiplic- 
ity covariance can be derived essentially independent of details of hadronization. The 
range of correlation in pseudo-rapidity emerges from the study without any dynamical 
"q ■ input, thus dispelling the notion that correlation length has any fundamental signifi- 

cance. An attempt is made to relate the results to the two-component structure seen 
in autocorrelation. 



(in 1 Introduction 

t> 

Correlations among particles produced at high energies have always been a subject of great 
^ \ interest, starting from the beginning in hadronic collisions [TJ, and more recently in nuclear 

collisions [2]. Over the years an abundant supply of experimental data have been accumu- 
lated, some of which have very high statistics, especially in the work done at relativistic 
heavy- ion collider (RHIC). By comparison, theoretical studies of the correlation phenom- 
ena in the bulk (excluding high-p^ jets) have been meager. Before the basic processes in 
hadron collisions mainly at low px are fully understood in a way that can be accepted as 
having a solid theoretical foundation, the profusion of nuclear data has inundated the subject 
that has no clear theoretical guidance, resulting in conflicting interpretations at times. The 
terms "long-range" and "short-range" correlations have been used, primarily in reference to 
what have been observed, rather than as properties of correlations understood at the level 
of parton interaction. 

STAR collaboration has produced extensive and detailed data on correlations in nuclear 
collisions [3l HJ EE El [7] as well as in pp collisions [SJ [9j [10] . The analyses have been done 
on autocorrelations in the angular difference variables = rji — r] 2 and 0a = 4>i ~ 02- 
The method used is to invert the scale-dependent (pr) fluctuations [TTJ [T2J, [13], yielding a 
rich structure in correlations not seen by any other method. Theoretical interpretation of 
the result is unfortunately grossly out of step with the growth of experimental information. 
Earlier models, such as the dual parton model [H], were adequate to treat the problem at 
a level commensurate with the coarseness of the data available at the time, but no recent 
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attempt has been made to calculate more accurately the basic correlations in pp collisions 
and to compare with the latest data. Soft production is still based on the idea of string 
fragmentation, the implementation of which is not much different from the original Lund 
model [15] . On nuclear collisions efforts have been made to incorporate string fusion [T6l fT7j. 
but no explicit calculation on correlation has been done. Very recently, the possibility of 
forward-backward correlation in the framework of color glass condensate has been advanced 
[18J , but without quantitative result yet to be compared with data. 

In view of the present theoretical status described above, a model-independent calcula- 
tion on any portion of the problem would lead to an improvement of that status. There 
is, however, no feasible formalism for calculating soft QCD processes involving multiparticle 
production at low px- Instead of starting from first principles, it is sensible to examine 
relevant hints from the data and restrict the scope to the issues most salient for correlations. 
It is furthermore reasonable to start with pp collisions before tackling the complexity asso- 
ciated with nuclear collisions. Careful analysis of the high-statistics data on the inclusive 
Pt distribution in pp collisions at y/s = 200 GeV has revealed a two-component picture 
of the particle production process essentially independent of theoretical models [19] . Infer- 
ences on the properties of parton scattering are made. The dominant soft component gives 
rise to a correlation behavior in rj/\ that is drastically different from the contribution from 
the minor, but significant, harder component, as revealed by the Porter- Trainor analysis 
of autocorrelation for pp collisions [HI El HO] . Although the autocorrelation measure is the 
Pearson's correlation coefficient that is the normalized version of the covariance of multiplic- 
ity densities, the t/a dependence is primarily the same as the forward-backward multiplicity 
covariance D"j b that has been determined without averaging over the rapidity sum rji + r^ 2 
[20]. 

We shall focus on just that piece of data on D^ b and relate it to the basic partonic 
scattering process that has back-to-back correlation. The connection between scattering 
angle and pseudo-rapidity is carefully done at the parton level to reveal the basic structure 
in the t?a dependence after hadronization is incorporated to relate the partonic properties 
to the hadronic observables. Since the processes are at low p?, there are inevitably some 
free parameters that are not fixed by our current understanding of soft physics. While the 
strength of correlation is adjustable to fit the data, there are no adjustable parameters to 
specify the range of correlation. Thus an understanding of the t/a dependence will be the 
main finding in this study that is essentially independent of models. Hadronization will be 
considered in the context of recombination, but details are not important in any essential 
way. A significant outcome of our study is the recognition that parton and hadron pseudo- 
rapidities can be sufficiently different that partons in one region may give rise to hadrons 
in another region, thereby cross-fertilizing forward and backward windows that are not too 
far apart. A correspondence of our result with the two components in the Porter- Trainor 
analysis will be attempted at the end of our study. 

2 Two-parton joint distribution 

We start with a brief discussion of the kinematical relationship between the parton and 
hadron variables. Since the forward-backward correlation (FBC) is measured in terms of 
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pseudo-rapidity without reference to the transverse momenta pr of the detected hadrons, it 
is clear that the correlation properties are dominated by the hadrons produced at low px, 
since high-p^ particles are severely suppressed. Since pQCD is not valid at low pt, there is 
no reliable theoretical formalism in which rigorous calculation can be made. Nevertheless, 
quarks and gluons are the constituents that interact, and those partons before and after 
scattering are the intermediary that bridges the initial state of the nucleons and the final 
state of multiparticles. We therefore seek some simple description of the parton state before 
hadronization in the hope that the observed FBC can be related to some basic properties of 
the partons. Since the hadronic px is not specified in the measurement of FBC, there is no 
point in emphasizing the partonic transverse momenta fcr, especially when their behavior 
is not calculable. Hence, our focus will be on the angular properties of the partons and 
hadrons. 

Let the parton and hadron momenta be denoted by k and p, respectively, and let their 
corresponding polar angles be 9 and G. If the angle between k and p is denoted by ip, then 
we have 

cos = cos 9 cos ip + sin 9 sin ip cos ip, (1) 

where ip is the solid angle between the plane containing k and p and the plane containing 
k and the z axis. The ip dependence in the problem is a property of the hadronization 
process. At low p? parton recombination is far more efficient, thus dominant, compared to 
fragmentation, which requires higher kj- partons. The dynamical mechanism need not be 
specified at this stage of geometrical consideration, except to mention that another parton 
beside the one at k, not necessarily collinear, has to be picked up to form the hadron at p. 
Let the hadronization cone be described by a Gaussian distribution in ip 



ip 2 

g(ip) = c/ exp(- — ), (2) 



where the normalization factor go is determined by 



2-7T pit/A 

d<p / dipsinip g{ip) = (N) . (3) 



o 

(N) is the average number of hadrons produced by a parton. The cone width o is expected 
to be small for pp collisions, and larger for AA collisions. For a much smaller than the range 
of angular integration in Eq. ([3]), we approximate sin?/' by ip so the integration can be readily 
carried out, yielding 



The differential form of Eq. ([3]) is 

dN 



d(pd cos ip 



gty). (5) 



At this stage of our discussion of hadronization it really does not matter whether hadrons 
are formed by recombination or fragmentation, since a hadronization cone is relevant in both 
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cases. We postpone further discussion about hadronization until the following section, since 
what is described above is sufficient to lead us from the hadrons to the partons. 

We now consider the two-parton distribution of back-to-back scattering at low kr in the 
CM system of the partons. Starting from the simplest possible form at this point, with 
allowance for complications to be included later, we write 

A S{z 1 + z 2 ), (6) 



dz\dz 



where z\ = cos^, assuming temporarily that the partonic rest frame is the same as the 
pp CM system. Since the two partons are exactly back-to-back, it is obviously a possible 
source of correlation at the partonic level. There can, of course, be other types of partonic 
interaction not describable in a simple form as in Eq. (J6j). However, it is not our problem here 
to provide a listing of all possible expressions of those interactions. Our aim is to investigate 
to what extent the observed FBC at the hadronic level can be traced back to Eq. ([6]). 

The 8{z\ + z 2 ) function is, of course, to constrain z 2 = —Z\ that is only valid if the CM 
systems of partons and hadrons coincide. In reality the proton has a wide distribution of 
partons at low x, so any pair of them from the two initial protons can have varying degree 
of mismatch between the two systems. Thus in the pp system we must allow 5{zi + z 2 ) to 
be broadened, but the way in which it is to be done is not calculable, since these are all 
low-Q 2 partons. We use a one-parameter description of the broadened distribution to replace 
5{z x + z 2 ): 



D m (zi + z 2 ) = — 



l-\{z x +z 2 f 



(7) 

which is forced to vanish at the kinematic limit \zi + z 2 \ = 2 and is normalized such that 



N, 



in 



J 2 dzD m {z) = 25(1/2, m + 1), 



B(a, b) being the Beta function. This distribution is roughly Gaussian shaped in the finite 
interval that Z\ + z 2 is allowed to vary, and can be broad when m is small. Since D is a 
constant, we know that in the limit m — > 0, D m {z\ + z 2 ) cannot give rise to any correlation 
between the two partons. Thus for the correlated part, we must subtract the uncorrelated 
part and introduce 

C m {z x + z 2 ) = D m {z x + z 2 ) - A), (9) 

which is what we shall use in place of 8{z\ + z 2 ) for the correlation to be calculated below. 
Let the pseudo-rapidity of a parton be denoted by (, i.e., 

C = -In tan 6/2, (10) 

from which can be derived 

z = cos 8 = tanh£, dz = cosh -2 £ d(. (11) 
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Using this in Eqs. (jSJ) and (J9]) yields 

dn c A 



d(id( 2 cosh d cosh ( 2 

where 



c m (G,C 2 ), (12) 



C TO (Ci,Ca) = ^[1 " ^(tanhCi + tanhC 2 ) 2 ] m - ~ (13) 

The subscript c in Eq. (j!2p is to emphasize that it is the correlated part of the two-parton 
distribution we are describing. 



3 Two-hadron joint distribution 

Having obtained an expression for the correlated two-parton distribution in d and (2, we 
now develop from it the correlated part of the two-hadron distribution. Equation ([T]) gives 
the relationship between the hadronic polar angle G in terms of the partonic polar angle 9, 
with t/> being the hadronization cone angle. Since soft hadrons with pr^l GeV/c dominate 
FBC, the participating partons that hadronize by recombination have transverse momenta 
hp in the range ^0.5 GeV/c. For a quark with k that hadronizes, it picks up an antiquark 
with momentum k! and forms a pion with momentum p = k + k'. The vectors k and k! need 
not be exactly collinear, but their angular difference cannot be too large, since the pion size 
is finite. In momentum space the uncertainty in the pion wave function is of the order of 
the pion mass, so the angular difference between k and k' is of order 0.14/0.5. The angle ip 
between k and p is about half that much. Thus the hadronization cone width a in Eq. (2) 
is roughly between 0.1 and 0.2. 

The hadronization process considered here implies that only one pion is produced per 
inclusive parton, since we are not discussing the total number of hadrons formed from a 
fixed pool of partons. Thus the normalization integral in Eq. ([3]) is (N) = 1. The precise 
value of (N) is unimportant in the situation where the parton density is not precisely known, 
exemplified by the undetermined parameter Aq in Eq. ([6]). We have not been thorough in 
our discussion about the recombination process, nor have we mentioned specifically the role 
that gluons play. The details of such processes, including the conversion of gluons to quarks 
before hadronization, have been studied extensively before [211 [22] . They are omitted here 
since such details would not contribute to the clarification of any crucial issues at hand and 
would only distract the flow of our main concern here. 

For a fixed parton momentum k(9) the probability of finding a hadron at angles (ip, ip) 
relative to it is given by Eq. (J5J). To guarantee that Eq. (pQ) is satisfied for the hadron at B, 
we define the hadron distribution per unit cos G as 

dN r 2n r n / 4 dN 

G(9, 6) = — = / dip / dip- — — <5(cos6 — cos 6* cos — sin # sin cos y?) (14) 

acos0 Jo Jo dipdyj 

such that 

J d cos eG(9,Q) = 1, (15) 
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where the range of integration corresponds to the angles in Eq. ([3]). Carrying out the 
integration over ip in Eq. (fl4j) we obtain 

G(6,B) = 2 J dip sinip g(ip)\ sin 2 6 sin 2 ip - (cos6 - cosflcos^) 2 ]- 1 / 2 . (16) 

In the realistic 3D geometry of the experiments the angular measure is in d cos dep. The 
azimuthal angle <p around the beam axis is integrated over 2n in both [TU] and [20], and is 
of no concern here. 

With the definition Z = cosG the hadron distribution in Z is related to the parton 
distribution in z as 

lz = J Tz (z ' } - ( } 

In terms of the hadron pseudo-rapidity rj, as in Eq. (TTTT) . 

Z = tanhr/, dZ = cosh -2 rj drj, (18) 

we have 



Defining 



f-/4^ GW0 ' w (19) 



G(C>7) = — ^-GWO.^W), (20) 
cosh 7/ 



we then have for the two-hadron distribution 

dN dn 



/(ITt 
dC *5c* G(Cl '" l)G(C2 '" 2) - (21) 

The hadronization function G((,r]) has no simple analytical form and is non-trivial. For 
small a, G((,r]) may be shaped roughly as a Gaussian, depending on (. Because of the 
non-linear relationship between polar angle and pseudo-rapidity, what is symmetric in ip as 
in Eq. (j2j) cannot be symmetric in 77 around (. In Fig. 1 we show some illustrative examples 
of G((,r]) for ( = 0.2 and a = 0.1,0.2 and 0.5. As a increases, the asymmetry in rj — ( 
develops in the wings, and becomes more significant for higher ( (not shown). The physics 
of hadronization is in the physical 3D momentum space involving ip, not in £ or rj. In AA 
collisions a can be large and G((, rj) can be wide in 77 — £. It is then unreliable to assume 
that a window in rj for the detected hadrons corresponds to a similar window in ( for the 
originating partons. 

Returning to the two-hadron distribution given in Eq. ff2TT) we note that the two partons 
at £1 and (2 hadronize independently, but because of the finite width in rj — £ of the G 
functions, it is possible that for any given window in rj the detected hadrons in any given 
event can originate from both partons separately at different d and £2- In other words, 
hadronization can "spill over" from positive to negative sides, and vice-versa, even if the 
forward and backward windows are located symmetrically on the two sides of rj = 0. 
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Putting Eq. (fT2|) in fl21|) we obtain for the correlated two-hadron distribution 



= / ^ co*' Cl c'osh' C 2 °^ ^ M ' »>• (22) 

Note that without the C m (£i,£2) function the integral is factorizable, resulting in no corre- 
lation. Thus the non-factorizable C m (d, (2) is the source of FBC. The inverse-square factor 
(cosh (1 cosh (2)^ suppresses the large ( contribution, so dN c /drjidr]2 appears to have no 
long-range correlation. However, that factor is the Jacobian of the transformation from z to 
( in Eq. (JTT1) and has no dynamical content. Thus the designation of such terms as short- 
or long-range correlation can be misleading. 



4 Forward-backward correlation 

To get FBC we integrate Eq. ( 1221) over 771, 772 with 771 in the forward window, and 772 in the 
backward window, defined to have widths £77 and spaced symmetrically apart from 77 = 
with 77A being the distance between the centers of the windows. Thus the correlated part of 
the forward and backward multiplicities is 

rv+ r~v- dN 
(N F N B ) C ( VA )= dm d V2 —±- 
Jr,- J-v+ dr]idri2 

= /<*&/ d( 2 , 2 /° , 2 , C m (Ci, C2)#iKCl, ^)H B (C,2, Va), (23) 
J J cosh Qi cosh (2 

where r\_ = (77^. — Srj)/2 and 77+ = (//a + ^)/2 , and 

h f (Cuva)= r dmGiCum), (24) 

Jr t - 

H b ((2,Va)= r~ d V2 G{( 2 ,r, 2 ). (25) 

It is clear in Eq. (|2"31 that hadronization is described entirely by the Hp and Hb functions. 
The hadronization of the two partons at (± and (2 can allow a forward (backward) parton into 
a backward (forward) window when 77A is small, since Hp and H b are not narrow functions. 

To illuminate the properties of hadronization, we show in Fig. 2 Hp{(, t?a) vs £ for fixed 
77A = 0.4 and £77 = 0.2, and for three values of a. All curves peak at ( = 7/a/2 = 0.2, which 
is the location of the window extending from 0.1 to 0.3, as indicated by the shaded interval. 
The solid curve shows that even for the narrow cone width of o = 0.1 partons significantly 
outside that window can contribute. It is even wider for the dashed curve. Although large a 
(such as 0.5) is not relevant to pp collisions, it shows that the value of ( of the contributing 
parton can be greater than 1, i.e., outside the detector coverage of STAR. 

The covariance of N f Nb that is measured is denoted by D"j b : 

D% = (N F N B ) - (N F ) (N B ) , (26) 

whose 7/a dependence in the data [20] is shown in two figures below. Identifying D"j b = 
(N f Nb) c with Eq. (1231) . we can fit the data by varying the relevant parameters in the 
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problem. Before discussing how that is done, it is important to note first that the decrease 
of D^ b with increasing 77 a is strongly affected by the (cosh d cosh C2) -2 factor in the integrand 
in Eq. ( J23l) . Changing that factor artificially to a weaker power diminishes the rate of decrease 
of (NfNb) c with 77A and will not fit the data. Since there is nothing to adjust in that factor 
to reflect the nature of interaction, the range of correlation that one may naively infer from 
the data has no dynamical meaning. 

There are essentially only two parameters to adjust in the problem: A Q and m. A 
describes the strength of parton interaction at low p F , and m specifies the spread of the 
parton rest frame relative to the pp CM system, thereby affecting the observed properties 
of correlation, a is between 0.1 and 0.2, and will not affect the final result sensitively. We 
adopt the following fitting strategy. First, we fix Aq — 1, a = 0.2, and vary m to see the 
corresponding dependence of (NfNb) c (m,7?A) on 7/a- That is shown in Fig. 3. Evidently, 
there is strong dependence on m. For narrower C m (d, (2) at higher m, the hadronic FBC is 
stronger, as it should for fixed A . However, the dependence on 7/a appears to be universal. 
We plot the ratios 

(N F N B ) c (m) 

RmiVA) ~ (N F N B ) e {m=l) (2?) 

in the inset and find that they are essentially independent of 7?a, except when m = 10 or 
larger. We do not expect D m to be narrow, so m should not be large. For smaller m, R m 
increases approximately linearly with m. A fit of the whole range of m studied results in 

R m = 0.0517 + 1.0335m - 0.1083m 2 + 0.005m 3 , (28) 

at r/A = 0.6. It is clear then to fit the data on D"j b the value of A must decrease with m as 
i?" 1 . Let us then define 

B m = A R m (r] A = 0.6), (29) 

which we expect to be a rather stable variable to use to fit Dj b for a wide range of m. 
Still holding a = 0.2, we obtain the results shown collectively by the solid line in Fig. 4 for 
m = 0.5, 1, 2. There is no dependence on m. The value of B m is 

B m = 14, m < 3. (30) 

The reproduction of the data in their tja dependence is very good, except for the highest point 
at 77a=1-6. We emphasize that we only vary B m to fit the normalization; the r/A dependence 
follows from Eq. ( |23l) without adjustment. For a = 0.1 we fix all other parameters already 
considered, and obtain the dashed line in Fig. 4 for m = 1. For other values of m there are 
some small variations because R m that is used is given in Eq. fT28|) . which was calculated 
for a = 0.2. But the difference is miniscule (< 5%). The dependence on a is evidently not 
strong, so we shall hereafter adhere to the value a = 0.2. 

To summarize our result so far, we see that A Q and m are strongly correlated in the fitting 
procedure. A C m (d, ( 2 ) is the source of hadronic correlation that replaces the A 5(zi + z 2 ) in 
Eq. <Q. The smaller m is, the larger must be A to compensate for the spread of C m (Ci, C2) 
in order to account for the observed magnitude of FBC. We have combined the two into 
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one parameter B m and succeeded in fitting the data. But the adjustment is only in the 
magnitude of the correlation. The important point is that the basic dependence of (N F N B ) c 
on 7?a is due to the factor (cosh (j cosh C2)~ 2 in Eq. (|23l) that is not adjustable and is the 
origin of the universality seen in Fig. 3. 

Let us, for definiteness, consider the case m — 1, for which B,\ — 1, so Aq — 14. Using 
this number in Eq. ([6]) and integrate it over —1 < z 2 ) < 1, we obtain 28 correlated pairs 
of partons. This is not a large number compared to the total number of pairs of partons 
(n(n — 1)), where (n) is roughly twice the total average charged event multiplicity (N) of 
hadrons, which is about 20. Thus only a small fraction (< 2%) of all parton pairs are 
correlated. 

Returning now to Eq. Q) that expresses the two-parton distribution assumed from the 
start, one may regard that as the s-wave contribution to the scattering process. It is of 
interest to ask whether the data on FBC can admit a p-wave contribution as well. To that 
end we extend Eq. §5§ to include a second term 

fin 

[A + A 1 (z 2 + z 2 2 )]5(z 1 + z 2 ), (31) 



dz\dz 2 



which is symmetric under z\ <-» ±z 2 . Since all quantities considered here and in the experi- 
ment are integrated over the 2n range of the azimuthal angle 0, the cross term between s and 
p waves vanishes by azimuthal symmetry. Applying Eq. ( !3~TT) to Eq. ( |23l) we replace the nu- 
merator by Aq + Ai(tanh ( 2 + tanh ( 2 ). The p-wave contribution increases with Q because of 
the factor (tanh (f + tanh £f) , so it may lift the upper end of the t/a distribution of (NfNb) c 
for a possibly better fit of the data. For values of Aq and A\ to be specified below, the 
contributions of the s- and p-wave components are shown by the dashed and dashed-dotted 
lines, respectively, in Fig. 5; the latter indeed increases with rj^. The solid line represents 
their sum, whose fit of the data is about comparable to that in Fig. 4. However, we have 
another motive for considering the additional p-wave contribution. 
To explain that, it is best to rewrite Eq. fl3T|) first as 

dTi 

[A s + A h {z 1 ,z 2 )]S(z 1 + z 2 ), (32) 



dz\dz2 
where 

A s = A -0AA 1 , (33) 

A h ( Zl ,z 2 ) = A 1 {0A + z 2 l +z 2 2 ). (34) 

The inclusion of 0.4Ai in Eq. ( 134|) is for the purpose of rendering the Ah contribution roughly 
constant in tja- That is shown by the dotted line in Fig. 5. The overall fit of the data is now 
achieved by choosing 

A = 13.2, A x = 0.2A = 2.64. (35) 

The dotted line is the contribution from £2) magnified by a factor of 2 for visual 

clarity, as is done for the dashed-dotted line. Without the multiplier, Ah is roughly A s /10 
at small 7/a; that is the constraint that leads to the determination of A\ relative to Aq in 



9 



Eq. (|35|) . Although the fit of the data on D"j b in Fig. 5 is not significantly improved, we have 
demonstrated that the data can be understood as a combination of two components, one of 
which is nearly flat in //a- 

The reason for doing the above decomposition is to establish a connection with the picture 
formed in the data analyzed in Ref. [19] that there exist a soft and a hard component; the 
former is larger and has a Gaussian-like decrease in rj^, while the latter is smaller and 
roughly constant in //a [HI El [10] . Our notations A s and are chosen to make a symbolic 
correspondence to those soft and hard components. The magnitude of the hard component 
depends on the charge multiplicity of the event class analyzed. The ratio of hard to soft 
components can be as much as 10% but can be much lower. The choice of A\ in Eq. ( 1351) 
and the solid line in Fig. 5 correspond to the maxium value that A\ can have. For weaker 
hard components, A\ is lower and the height of the dotted line is also lower accordingly. The 
quality of the overall fit of the data on D 2 j b in Fig. 5 is largely unaffected. 

The soft and hard components in Ref. [19] are determined by detailed analysis of the pt 
behavior of the produced hadrons, the hard having (p T ) only around 1.2 GeV/c, but still 
higher than the soft component. Our study here bypasses the issue about the transverse 
momenta of the partons and hadrons, which are not tractable at such low p^, but focuses on 
the angular variables that are more intimately related to pseudo-rapidity. We find that the 
properties found in Refs. [HI [9j [10j [19] can be accommodated by a combination of s- and p- 
wave components to give the A s and Ah components. Since we have had no need to question 
the kx aspect of the partons in the present treatment, we have no basis to judge what is soft 
or hard. However, it does seem reasonable that the soft component consists entirely of the 
s-wave, while the hard component includes also the p-wave part of the scattering. 

5 Conclusion 

We have considered hadronic FBC by relating it to the simplest form of partonic interaction 
that has back-to-back correlation. It is found that it gives a fairly good description of 
the dependence of FBC on window separation rj^ without any adjustment of the range 
of correlation. The magnitude of the correlation function is fitted by varying the number 
of correlated pairs of partons, which turns out to be less than 2% of all possible pairs 
that hadronize. The so-called "range of correlation" in rj is mainly a consequence of the 
transformation from the polar angles, and has no meaning at the parton level. 

Since low-pr hadrons dominate any hadronic measure that does not restrict the pt range 
of coverage, it is necessary to consider \ow-kr partons for which available theoretical tools are 
deficient. Leaving open the questions about kx and pr, we have focused on the relationship 
between polar angles and pseudo-rapidities and found a meaningful way to separate different 
issues that include the mismatch between the partonic and hadronic CM systems in the initial 
state, and the hadronization of partons in the final state. Uncertain properties of those two 
specific issues have been investigated, and various possibilities have been considered before 
arriving at the final result that is mainly insensitive to those properties. 

Since at pt < 1.5 GeV/c recombination may be more likely than fragmentation as the 
dominant mechanism of hadronization, we have considered the effect of non-vanishing width 
of the hadronization cone. We find that there is significant cross feeding of partons from 
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one region to hadrons in neighboring regions. To identify the partonic window with the 
hadronic window is an assumption that becomes even more unreliable in nuclear collisions, 
where attempts have been made to infer the nature of the bulk medium from the observed 
charge fluctuation in restricted windows. 

We have attempted to make contact with the results of Porter- Trainor analysis that 
shows by autocorrelation the existence of very different t/a behaviors for the soft and hard 
components at low px (< 1.5 GeV/c) [81 [9J, [10] • We can identify the (soft) component having 
strong 7]/± dependence with our s-wave component and the weaker (hard) component having 
roughly no dependence on r/^ with our part that includes the p-wave contribution. Our 
major finding is that the strong r)&, dependence does not imply short-range correlation. It 
is possible that the ^-independent component may be related to multiplicity fluctuation in 
some way by treating the "minijets" in the hard component at the parton level. 

The parton-based approach adopted in this work is clearly different from the string model 
that emphasizes the role of the valence quarks in the formation of strings and the mechanism 
of string fragmentation for soft production of particles. Our approach is closer to the original 
parton model, where low-momenta soft partons are a part of the initial state of the incident 
proton [23] . 

With the elucidation that we have achieved for the pp collision problem, it is natural to 
ask what light it sheds on the AA collision problem. To the extent that we have investigated 
the subject we have found that the nuclear problem is vastly more complicated, since there 
are many contributing factors that can influence FBC. What we have uncovered in the 
pp problem is likely to be overwhelmed by fluctuations related to the particle production 
processes outside the realm of partonic correlation in pp collisions. It is therefore mostly a 
separate problem only a part of which is connected with the topic of study here. 
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Figure 2: The distributions of H F ((,r]^) for a forward window at = 0.2 and window 

size 5r) = 0.2 indicated by the shaded interval. 
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Figure 3: Examples of the FBC function in t]^ for fixed A — 1 and a = 0.2, and for six 
values of m that characterizes the width of C m (Ci, C2)- Ratios of those functions relative to 
m — 1 are shown in the inset. 
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Figure 4: Forward-backward multiplicity covariance versus t/a- Data are from Ref. [20]. The 
solid line shows the coalescence of three cases for a = 0.2 and m = 0.5, 1, 2. The dashed line 
is for a = 0.1 and m = 1; other cases for m = 0.5 and 2 are very nearly the same. 



16 




0.5 1 1.5 2 

Figure 5: Contributions to D"j h from the s-wave component in dashed line and p-wave com- 
ponent in dashed-dotted line, raised by a factor of 2 for visual clarity. The solid line is their 
sum. The dotted line represents the (hard) component specified by Ah in Eq. (I34p . also 
magnified by 2. The data are from [20] . 
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